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Abstract. This paper considers the problem of minimizing the ordered weighted average (or ordered 
median) function of finitely many rational functions over compact semi-algebraic sets. Ordered weighted 
averages of rational functions are not, in general, neither rational functions nor the supremum of rational 
functions so that current results available for the minimization of rational functions cannot be applied to 
handle these problems. We prove that the problem can be transformed into a new problem embedded 
in a higher dimension space where it admits a convenient representation. This reformulation admits a 
hierarchy of SDP relaxations that approximates, up to any degree of accuracy, the optimal value of those 
problems. We apply this general framework to a broad family of continuous location problems showing 
that some difficult problems (convex and non-convex) that up to date could only be solved on the plane 
and with Euclidean distance, can be reasonably solved with different ^p-norms and in any finite dimension 
space. We illustrate this methodology with some extensive computational results on location problems in 
the plane and the 3-dimension space. 

1. Introduction 

Weighted Averaging (OWA) or Ordered Median Function (OMF) operators provide a parameterized 
class of mean type aggregation operators (see [25] [44] and the references therein for further details) . Many 
notable mean operators such as the max, arithmetic average, median, k-centrum, range and min, are mem- 
bers of this class. They have been widely used in location theory and computational intelligence because 
of their ability to represent flexible models of modern logistics and linguistically expressed aggregation in- 
structions in artificial intelligence ([35] and [321 EH EH [42 , 43, 44 ). Weighted averages (or ordered median) 
of rational functions are not, in general, neither rational functions nor the supremum of rational functions 
so that current results available for the minimization of rational functions are not applicable. In spite 
of its intrinsic interest, as far as we know, a common approach for solving this family of problems is not 
available. Nevertheless, one can find in the literature different methods for solving particular instances of 
problems within this family, see e.g. [5J El M 123 US HH M 1301 133 131 E3 • The first goal of this paper 
is to develop a unified tool for solving this class of optimization problems. In this line, we prove that the 
general problem can be transformed into a new problem embedded in a higher dimension space where it 
admits a convenient representation that allows to arbitrarily approximate or to solve it as a minimization 
problem over an adequate closed semi-algebraic set. Hence, our approach goes beyond a trivial adaptation 
of current theory. 

Regarding the applications, it is commonly agreed that ordered median location problems are among the 
most important applications of OWA operators. Continuous location has achieved an important degree of 
maturity. Witnesses of it are the large number of papers and research books published within this field. In 
addition, this development has been also recognized by the mathematical community since the AMS code 
90B85 is reserved for this area of research. Continuous location problems appear very often in economic 
models of distribution or logistics, in statistics when one tries to find an estimator from a data set or in 
pure optimization problems where one looks for the optimizer of a certain function. For a comprehensive 
overview the reader is referred to [4j or |25j . Despite the fact that many continuous location problems 
rely heavily on a common framework, specific solution approaches have been developed for each of the 
typical objective functions in location theory (see for instance [J]). To overcome this inflexibility and to 
work towards a unified approach to location theory the so called Ordered Median Problem (OMP) was 
developed (see |25j and references therein). Ordered median problems represent as special cases nearly 
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all classical objective functions in location theory, including the Median, CentDian, center and k-centra. 
More precisely, the 1-facility ordered median problem in the plane can be formulated as follows: A vector 
of weights (Ai, . . . , A n ) is given. The problem is to find a location for a facility that minimizes the weighted 
sum of distances where the distance to the closest point to the facility is multiplied by the weight A„ , the 
distance to the second closest, by A n _i, and so on. The distance to the farthest point is multiplied by 
Ai. Many location problems can be formulated as the ordered 1-median problem by selecting appropriate 
weights. For example, the vector for which all Aj = 1 is the unweighted 1-median problem, the problem 
where A„ = 1 and all others are equal to zero is the 1-center problem, the problem where Ai = . . . = A& = 1 
and all others are equal to zero is the /c-centrum. Minimizing the range of distances is achieved by Ai = 1, 
A n = — 1 and all others are zero. Despite its full generality, the main drawback of this framework is the 
difficulty of solving the problems with a unified tool. There have been some successful approaches that are 
now available whenever the framework space is either discrete (see [H [3U] ) or a network (see [TT] , [T^] 
or [H]). Nevertheless, the continuous case has been, so far, only partially covered. There have been some 
attempts to overcome this drawback and there are nowadays some available methodologies to tackle these 
problems, at least in the plane and with Euclidean norm. In Drezner [3j and Drezner and Nickel [6J the 
authors present two different approaches. The first one uses a continuous branch and bound method based 
on triangulations (BTST) and the second one on a D-C decomposition for the objective function that allow 
solving the problems on the plane. More recently, Rodriguez-Chia et al. [51] also address the particular 
case of the /c-centrum problem and using geometric arguments develop a better algorithm applicable only 
for that problem on the plane and Euclidean distances. 

Quoting the conclusions of the authors of [5]: "All our experiments were conducted for Euclidean 
distances. As future research we suggest to test these algorithms on problems (even the same problems) 
based on other distance measures. (...) Solving k-dimensional problems by a similar approach requires the 
construction of k-dimensional Voronoi diagrams which is extremely complicated. " 

Therefore, the challenge is to design a common approach also to solve the above mentioned family of 
location problems, for different distances and in any finite dimension. This is essentially the second goal 
of this paper. In our way, we have addressed the more general problem that consists of the minimization 
of the OWA operator of a finite number of rational functions over closed semialgebraic sets that is the 
first goal of this paper. Thus, our second goal is to solve a general class of continuous location problems 
using the general approach mentioned above for the minimization of OWA rational functions and to 
show the powerfulness of this methodology. Of course, we know that the problem in its full generality is 
NP — hard since it includes general instances of convex minimization. Therefore, we cannot expect to 
obtain polynomial algorithms for this class of problems. Rather, we will apply a new methodology first 
proposed by Lasserre |16j , that provides a hierarchy of semidefinite problems that converge to the optimal 
solution of the original problem, with the property that each auxiliary problem in the process can be solved 
in polynomial time. 

The paper is organized in 5 sections. The first one is our introduction. In the second section and 
for the sake of completeness, we recall some general results on the Theory of Moments and Semidefinite 
Programming (SDP) that will be useful in the rest of the paper. Section [3] considers what we call the 
MOMRF problem which consists of minimizing the ordered median function of finitely many rational 
functions over a compact basic semi-algebraic set. In the spirit of the moment approach developed in 
Lasserre [16j HH] f° r polynomial optimization and later adapted by Jibetean and De Klerk [10] , we define 
a hierarchy of semidefinite relaxations (in short SDP relaxations). Each SDP relaxation is a semidefinite 
program which, up to arbitrary (but fixed) precision, can be solved in polynomial time and the monotone 
sequence of optimal values associated with the hierarchy converges to the optimal value of MOMRF. 
Sometimes the convergence is finite and a sufficient condition permits to detect whether a certain relaxation 
in the hierarchy is exact (i.e. provides the optimal value), and to extract optimal solutions (theoretical 
bounds on the relaxation order for the exact results can be found in [35j[36]). Section [4] considers a general 
family of location problems that is built from the problem MOMRF but which does not actually fits under 
the same formulation because the objective functions are not quotients of polynomials. Nevertheless, we 
prove that under a certain reformulation one can define another hierarchy of SDP that fulfils convergence 
properties 'a la Lasserre'. This approach is applicable to location problems with any £ p -norm (p 6 Q) 
and in any finite dimension space. We exploit the special structure of these problems to find a block 
diagonal reformulation that reduces the sizes of the SDP relaxations and allows to solve larger instances. 
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Our computational tests are presented in Section 5. We analyze five families of problems, namely, Weber, 
center, fc-centrum, trimmed-mean and range. There we show that convergence is rather fast and very high 
accuracy is achieved in all cases, even with the first feasible relaxation. (We observe that for location 
problems with Euclidean distances that relaxation order is r = 2.) The paper ends with some conclusions 
and an outlook for further research. 

2. Preliminaries 

In this section we recall the main definitions and results on the moment problem and semidefinitc 
programming that will be useful for the development through this paper. We use standard notation in the 
field (see e.g. gU]). 

We denote by R[x] the ring of real polynomials in the variables x = (xi, . . . , x n ), and by K[x]d C R[x] 
the space of polynomials of degree at most d £ N (here N denotes the set of nonnegative integers). We 
also denote by B = {x a : a £ N™} a canonical basis of monomials for R[x], where x a = x" 1 ■ ■ •x" n , for 
any a£N" 

For any sequence indexed in the canonical monomial basis B, y = (y a )a£fi n C R, let L y : R[x] — > R be 
the linear functional defined, for any / = 2 aGN „ fa x " € M[x], as L y (/) := 2 q(E n» fa Vol- 

The moment matrix M ( j(y) of order d associated with y, has its rows and columns indexed by (x a ) and 
M d (y)(a,y9) := L y (x a+l3 ) = y a+ p, for |a|, |/?| < d. Note that the moment matrix is ( Tl + d ) x and 
that there are (™~^ 2d ) y a variables. 

For g £ R[x] (= 57gN™x 7 ), the localizing matrix M<j(<7, y) of order d associated with y and g, has its 
rows and columns indexed by (x a ) and Md(g,y)(ct, f3) := L y (x Q+ ' 3 g(x)) — J2 7 Q-yV-y+a+Pi for \a\, \/3\ < d. 

Definition 1. Let y = (y Q ) C R be a sequence indexed in the canonical monomial basis B. We say that y 
has a representing measure supported on a set K C R n if there is some finite Borel measure fi on K such 
that 

y a = x a dfi(x), for all a £ R". 

The main assumption that is needed to impose when one wants to assure the convergence of the SDP 
relaxations for solving polynomial optimization problems (see for instance [THJ HO]) was introduced by 
Putinar [33] and it is stated as follows. 

Putinar's Property. Let {gi, . . . ,gi} C R[x] and K := {x £ R d : gj(x) > 0, : j = 1, . . . ,£} a basic closed 
semialgebraic set. Then, K satisfies Putinar's property if there exists b £ such that: 

(1) {x : u{x) > 0} C R™ is compact, and 

(2) u = Co + 2j=i a j 9j> f or some f i, . . . , <7\ £ Y\x\. (This expression is usually called a Putinar's 
representation of u overK). 

Being S[x] C R[x] the subset of polynomials that are sums of squares. 

Note that Putinar's property is equivalent to impose that the quadratic polynomial M — JDILi x l has 
a Putinar's representation over K. 

We observe that Putinar's property implies compactness of K. It is easy to see that Putinar's property 
holds if either {x : gj(x) > 0} is compact for some j, or all <?j are affine and K is compact. Furthermore, 
Putinar's property is not restrictive at all, since any semialgebraic set K for which is known that Y^7=i x 1 — 
M holds for some M > and for all x £ K, K = K U {gi+i(x) := M - Yh=i x i > °l verifies Putinar's 
property. 

The importance of Putinar's property stems from the following result: 

Theorem 2 (Putinar [33J). Let {g u . . . ,gi} C R[x] and K := {x £ R d : gj(x) > 0, : j = 1, . . . , 1} satisfying 
Putinar's property. Then: 

(1) Any f £ R[x] which is strictly positive on K has a Putinar's representation over K. 

(2) y = (y a ) has a representing measure on K if and only «/M<j(y) ^ 0, and Md(<7j,y) ^ 0, for all 
j = 1, . . . ,1 and d£N. 

(Here, the symbol y stands for semidefinite positive matrix.) 
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The following result that appears in [TO] and [Kj will be also important for the development in the next 
sections. 

Lemma 3. Let K C M. d be compact and let p, q be continuous with q > on K. Let A4(K.) be the set of 
finite Borel measures on K and let 'P(K) C A4(K) be its subset of probability measures on K. Then 



J K pd/i 



min { / 



pdip : / qd<^ = 1} 



K 



mm 

M e-P(K) 



— d/i 
K <7 



aim —— — . 
c6K q[x) 



3. Minimizing the ordered weighted average of finitely many rational functions 
Let K C l d be a basic semi-algebraic set defined as 

K := {xeW i :g j (x)>0, j = l,...J} 

for g u ...,gt€ R[x}. 

Let us introduce the function OM(x) = Y^k=i ^k{ x )f(k)( x )i f° r some rational functions (fj) C M.[x], 
being f k = p k /q k rational functions with p k ,q k e R[x], A fe (x) G and f( k) {x) e • • • , fm(x)} 

such that /(i)(x) > /(2)(#) > *•• > f(m)( x ) f° r ^ G K n . We assume that K satisfies Putinar's property 
and that % > on K, for every k = 1, . . . , m. 

Consider the following problem: 

(OMRP A ) p A := min{OM(x) : x e K }, 



1. 



Associated with the above problem we introduce an auxiliary problem. For each i — 1, . . . , m, j 
. . . , m consider the decision variables Wa that model for each x € K 



Now, we consider the problem: 
(OMRP A ) 



w, 



1 if fi(x) = f(j)(x), 
otherwise. 



mm 

x.w 



^ ^j{x)^2fi{x)w i:j 



3=1 



(1) 



S.t. 



E* 

3=1 
m 

E' 
<4- 



= 1, for i = 1, 



)ij = 1, for J = 1, . 
Wy = 0, for i,j = 



i TO, 



(2) 



^Wijfiix) > X/ 



i=l 



2 = 1 



J = 1, ■••,to, 



(3) 
(4) 



Wij < TO, 



EE' 

»=i j=i 

lOjj € R, for i, j = 1, . . 



, m, x G K. 



The first set of constraints ensures that for each x, fi{x) is sorted in a unique position. The second 
set ensures that the j th position is only assigned to one rational function. The next constraints are 
added to assure that Wij € {0,1}. The fourth one states that /(i)(x) > ••• > f( m )(x). The last set of 
constraints ensures the satisfaction of Putinar's property of the new feasible region. (Note that this last 
set of constraints are redundant but it is convenient to add them for a better description of the feasible 
set.) 

These two problems, (OMRP A ) and (OMRP A | satisfy the following relationship. 



Theorem 4. Let x be a feasible solution of (OMRP A | then there exists a solution (x,w) for (OMRP A 



such that their objective valu es are eq ual. Conversely, if (x,w) is a feasible solution for (OMRP A | then 
there exists a solution x for (OMRP A I having the same objective value. Ln particular g\ = g\. 



MINIMIZING ORDERED WEIGHTED AVERAGING OF RATIONAL FUNCTIONS 



5 



Proof. Let a; be a feasible solution of (OMRP" ). Then, it clearly satisfies that x G K. In addition, let a 
be the permutation of (1, . . . , to) such that f a (i)( x ) > fa(2)i x ) > ■ ■ • > fa(m)( x )- Take, 

1 iii = a(j), 
otherwise. 

Clearly, (x,w) satisfy the constraints in (lpl- Indeed, for any i 53j=i ^ij = ^io- 1 ^) = 1- Analogously, 
for any j, Y^hLi ™ij = ^<r(J)J = !• By its own definition, u; only takes 0, 1 values and thus, wfj — iy = 
for all i,j and YITj ®fj — m - Finally, to prove that (a;, w) satisfies (p|, we observe, w.l.o.g., that for any j 



i=l 



there exist i* and i such that er(j') = i* and er(j + 1) = i. Hence,: 

m m 

J2 w vM x ) = w^jf^ix) > w ij+1 f a{j+1) (x) = ^w lJ+1 f l {x). 

i=l 

Moreover, 

m m 
3=1 «=1 



Conversely, if (x, iu) is a feasible solution of ( OMRPa I then, clearly x is feasible of (|OMRP^|) and by 



the above, OM x (x) = J2T=i ^j( x ) Yh=\ M x ) 1 



a 



Then, we observe that fi = Villi f° r each i — l,...,m. Therefore, the constraint X)i=i w ijfi( x ) ^ 
Y^JLi w ij+ifi( x ) can be written as a polynomial constraint as 

m m m m 

^tUyPi(x) JJgfc(ar) < ^ Wy+i JJ g fc (a;) j = l,...,m. 

i—l k^i i—l k^i 



Let us denote by K the basic closed semi-algebraic set that defines the feasible region of (OMRPa 



Lemma 5. I/K C M m satisfies Putinar's property then K C K" +m satisfies Putinar's property. 

Proof. Since K satisfies Putinar's property, the quadratic polynomial x i-> := M — ||x||| can be 

written as u(x) = o~o(x) + 'Y^j=i a j( x )9j{ x ) for some s.o.s. polynomials (o-j) C T,[x\. Next, consider the 
polynomial 

m rn 

(x,w) i ^ r(x,u>) = M + m- \\x\\\ "XlX^r 

i=l 3=1 

<m : r(x, z) > 0} C M n + m is compact and moreover, r can be written 



Obviously, its level set {(x, w) £ 
in the form 



g(x 1 w) defining K 
. . 



r(x,w) = a (x) + ^2aj(x)gj(x) + lx (to- ^^w^-), 

3=1 *=1 3 = 1 

for appropriate s.o.s. polynomials (c^-) C £[a;,u>]. Therefore K satisfies Putinar's property, the desired 
result. 

□ 



Now, we observe that the objective function of (OMRPa ) can be written as a quotient of polynomials 
in E[x,w]. Indeed, take 



(5) 

Then, 
(6) 



p\(x,w) = }^Xj (x) Wjj pi (x) If gfc (a;) and q\(x,w) = q k (i 

j — 1 i—l k^i k — 1 



3=1 



i=l 



_ Px( x ,w) 
q\{x,w) 
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Then, we can transform Problem ( OMRPa I in an infinite dimension linear program on the space of 
Borel measures defined on K. 

Proposition 6. Let K C R" +m be the closed basic semi-algebraic set defined by the constraints \1 
Consider the infinite- dimensional optimization problem 

V\: p\ = min I / p\d\i : I q x dp = 1, fj, € M(K) } , 
x ' w UK Jk ) 

being p\, q\ 6 K.[aj, w] as defined above. Then p\ = p\. 



Proof. It follows by applying Lemma [3] to the reformulation of ( OMRPa I with the objective function 



written using p\ and q\ in ^ . □ 

The reader may note the great generality of this class of problems. Depending on the choice of the 
polynomial weights A we get different classes of problems. Among then, we emphasize the important 
instances given by: 

(1) A = (1, 0, . . . , 0, 0) which corresponds to minimize the maximum of a finite number of rational 
functions, 

(2) A = (1, 1,0, ... ,0) which corresponds to minimize the sum of the fc-largest rational functions 
(fc-centrum) 

(3) A = (0, { k }), 0, 1, . . . , 1, 0, v. 2 ?, 0) which models the minimization of the (fei, /^-trimmed mean of m 
rational functions,... 

(4) A = (1, a, . . . , a) which corresponds to the a-centdian, i.e. minimizing the convex combination of 
the sum and the maximum of the set of rational functions. 

(5) A = (1, . . . , —1) which corresponds to minimize the range of a set of rational functions. 



Remark 7. Problem 



OMRPK 



can be easily extended to deal with the minimization of the ordered median 



function of a finite number of other ordered median of rational functions. The reader may observe that this 



can be done by performing a similar transformation to the one in (OMRPa) and thus lifting the original 
problem into a higher dimension space. 

3.1. Some remarkable special cases. 

The above general analysis extends the general theory of Lasserre to the case of ordered weighted 
averages of rational functions. Notice that this approach goes beyond a trivial adaptation of that theory 
since ordered weighted averages of rational functions are not, in general, neither rational functions nor 
the supremum of rational functions so that current results cannot be applied to handle these problems. 
However, one can transform the problem into a new problem embedded in a higher dimension space 
where it admits a representation that can be cast in the minimization of another rational function in a 
convenient closed semi- algebraic set. Needless to say that the number of indeterminates increases with 
respect to the original one. This may become a problem in particular implementations due to the current 
state of semidefinite solvers. 

In some important particular cases that have been extensively been considered in the field of Operations 
Research the above approach can be further simplified as we will show in the following. One of this cases, 
the minimization of the maximum of finitely many rational functions, has been already analyzed by Laraki 
and Lasserre [15] . We will show that the approach in [15] is also a particular case of the analysis that we 
present in the following. 

For the rest of this subsection we will restrict ourselves, for the sake of readability, to the case of scalar 
(real) lambda weights. We will begin with the case of A = (1, y?l, 1, . . . , 0), for 1 < k < m. Note that 
for the case k = 1 we will recover the case analyzed in [15j ■ the case k = m is trivial since it reduces to 
minimize the overall sum and the remaining cases are not yet known. 

We are interested in finding the minimum of the sum of the fc-largest values {fi(x), . . . , f m (x)} for all 
x € K, being a closed basic semi-algebraic set. In other words, for any k, k = 1, . . . , m — 1, we wish to 
solve the problem: 

k 
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We observe that for a given x, we have: 

k mm 

s k{x) = ^2f(j)( x ) = ma - x {^2 v jfj(x) ■ ^Z v j =k,0<Vj < 1, Vj}. 
i=i j'=i i=i 

Therefore, by duality in linear programming: 

m 

Sk(x) = minjfct + ■ : t + Tj > fj(x), rj > 0, Vj}. 
i=i 

Finally, we consider the problem: 

m 

g := min kt + rj 

j=l 

(kC) 8.t.t + rj>fj(x), j = l,...,m 

rj>0, j = l,...,m, 

ieK. 



Let us denote by K the basic closed semi- algebraic set that defines the feasible region of (kC). 

Lemma 8. IfK.CZ K™ satisfies Putinar's property then K C satisfies Putinar's property. More- 

over g = g. 

Proof. Since we have assumed K to be compact, for any j = 1, . . . , m, there exist LBj, UBj such that for 
any x £ K, 

LBj < fj(x) < UB r 

Let us denote LB = min J= i.. m LBj and UB = max J= i.. m UBj. Consider an arbitrary k, I < k < m — 1 
and an arbitrary (but fixed) x € K. Without loss of generality, assume that f m (x) > ... > fi{x). We 
define the function 

in 

g(t) := mm{kt + y"V j : t + r j > fj(x), r 3 > 0, V j = 1, ..,m}. 

3 = 1 

Clearly, g is piecewise linear and convex; and it attains its minimum on any point of the interval Ik = 
(fk+i(x), fk(x)]. Indeed, observe that for any t E Ik, the slope of g (i.e. its derivative with respect to t) is 
null since: 

k k 

g(t) = kt + J^ifM ~t) = J2 fi& = S *W- 

J = l 3=1 

^From the above, we observe that 

m 

g — min Sk(x) = min min{<7(i) : kt + rj :t + Vj > fj(x), Tj > 0, V j = 1, .., in} = g. 



It remains to prove that K, the feasible region of problem (kC), satisfies Putinar's condition. First, 
we observe from the argument above that in order to obtain the minimum value of the function g, for 
any k = 1, ..,m — 1 and any x 6 K, we only need to consider the range t £ (/( m )(x), f(u(x)]. Hence, the 
overall range for t can be restricted to LB < t < UB. On the other hand, for any x £ K, the constraints 
< r j < fj{x) — t set the range of the variable rj. Hence 

< Tj < UB 3 - LB, V j = 1, . . . , m. 

Including the constraints, LB < t < UB, < rj < UBj — LB, V j = 1, . . . , to, in the definition of K does 
not change the value of g and makes the feasible set compact. Thus, satisfying Putinar's condition. □ 
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This approach extends also to the more general case of non-increasing monotone lambda-weights, i.e. 
Ai > A2 > ... > A m > A„ l+ i := (Note that we define an artificial X m +i to be equal to 0). In this case 
the problem to be solved is: 



Q\ := mmMOM x (x) :=J2 X if(j)( x )- 

xE 3=1 

We observe that for a fixed x € K, we can write the objective function as: 

m 

3=1 

Then, we introduce the problem 

m 

(7) g x := min ^(A fc - \ k +i)S k (x) 

k=i 

tk + r kj > fj(x), j,k = l,...,m, 
rk 3 > 0, j,k = l,...,m, 
xeK. 



Let us denote by K the basic closed semi-algebraic set that defines the feasible region of the Problem 
([7]). Now, based in the previous lemma, it is straightforward to check the following result. 

Lemma 9. If K C K" satisfies Putinar's property then K C K n+m +m satisfies Putinar's property. 
Moreover gx = g\. 

Another class of problems that can also be analyzed giving rise to a more compact formulation that the 



one in the general approach ( OMRP>, I is the trimmed mean problem. A trimmed mean objective appears 

ki k-2 

for A = (677??7o, 1, . . . , 1, (WTO). 

This family of problems has attracted a lot of attention in last times in the field of location analysis 
because of its connections to robust solution concepts. Its rationale rests on the trimmed mean concepts in 
statistics where the extreme observations (outliers) are removed to compute the central estimates (mean) 
of a sample. Thus, we are looking for a point x* that minimizes the sum of the central functions, once we 
have excluded the &2 smallest and the hi largest. Formally, the problem is: 



n— &2 



i=kx + l 



Now, we observe that Yl7=k'+i f(i)( x ) ~ S n -k 2 ( x ) ~ Ski i x )- Therefore, using the above transformation 
we have: 

rn m 

SkA x ) = = max{y%j/j(>) - ^Vj = h,0 <Vj < 1, Vj}, 

3=1 3=1 
m 

S n -k 2 (x) = min{(n - k 2 )t + Y2 r i :t + r J - fj( x )> r j ^ °j V j}- 

3=1 
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Thus, using both reformulations the trim-mean problem results in: 
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m m 



(kTr) 



g := min (n - k 2 )t + ^2 r j - v ofj( x ) 

s.t. Vj = ki, 

t + rj>fj(x), j = l,...,m, 
r 3 >0, j = l,...,m, 
V j( v j- 1 )= Q , 3 = 1, 

ieK. 



Let us denote by K the basic closed semi- algebraic set that defines the feasible region of (kTr 



Lemma 10. If K c R™ satisfies Putinar's property then K C M™+ m + 1 satisfies Putinar's property. 
Moreover g = g. 



Remark 11. We observe that the special formulations for k-centrum (kCl and trim-mean (kTr 
specially suitable for handling these two classes of problems 

problem reduces to a ki-centrum, variables Vj are not needed and formulation (kTr 



First of all, we note that if k\ = the 
simplifies exactly to 



(kTr 



Second, we point out that both formulations take advantage of the special structure of the considered 

applied to these problems. 



problems and thus they are simpler than the general formulation ( OMRP 



Actually, the number of variables in (kC), for solving the k-centrum problem (resp. (kTr 
trim-mean problem), is m + n 



for solving the 



using (OMRP A 



1 (resp. 2m + d + 1 ) while the number of variables for the same problem 
This reduction is remarkable due to the current status of SDP solvers which 
are not at a professional level. In spite of that, those problems, where no special structure is known or 
it cannot be exploited, can also be tackled using the general formulation ( OMRP> ) at the price of using 
larger number of variables. 



3.2. A convergence result of semidefinite relaxations 'a la Lasserre'. 

We are now in position to define the hierarchy of semidefinite relaxations for solving the MOMRF 
problem. Let y = (y a ) be a real sequence indexed in the monomial basis (x^w 1 ) of W[x, w] (with a = 
(/3,7) e N™ x N™ 2 ). Let p x (x,w) and q x (x,w) be defined as in 

Let ho(x,w) := p\(x,w), and denote := [~(deg<7j-)/2~|, Vj := [~(deg/ij)/2] and Vj := f(deg^)/2] 
where {g\, . . . , gi} are the polynomial constraints that define K and {hi, . . . , h m } and {h[, . . . , h' m } are, 
respectively, the polynomial constraints ^ and ^ in K \ K, respectively. 

Let us denote by 1(0) = {l,...,n} and I(j) = {(j, k)}k=i,..., m , for all j = l,...,m. With x(I(0)), 
w(I(J)) we refer, respectively, to the monomials x, w indexed only by subsets of elements in the sets 
1(0) and I(j), respectively. Then, for with k = 1, ...,£, let M r (y, 1(0)) (respectively M r (gky, 1(0))) 
be the moment (resp. localizing) submatrix obtained from M. r (y) (resp. M r (gky)) retaining only those 
rows and columns indexed in the canonical basis of R[a;(/(0))] (resp. R[x(i~(0))]). Analogously, for hj 
and hp j = 1, . . . , m, as defined in (|2j) and (|3j), respectively, let M r (y, 1(0) U I(j) U I(j + 1)) (respectively 
M r (hjy,I(0)Ul(j)Ul(j + l)), M r (hjV,I(Q)Ul{j)Ul{j + l)) ) be the moment (resp. localizing) submatrix 
obtained from M r (y) (resp. WL r (hjy), M r (/i^y)) retaining only those rows and columns indexed in the 
canonical basis of R[x(I(0)) U w(I(j)) U w(I(j + 1))] (resp. R[x(J(0)) U w(I(j)) U w(I(j + 1))]). 
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For r > max{r , isq} where r 
ing hierarchy of semidefinite programs: 

min L y (p A ) 
y 

s.t 



(Qr) 



max f fc , v a 

k=l,...,£ 



max{ max v.j , max z/ } , we introduce the follow- 



j=i 



M r (y,/(0)) 




y 


0, 








U r ^ k (g k y,I(0)) 




y 


0, 


k = 


1,.. 


.,£, 


M r (y,/(0)U7(j)U/(i + l)) 




y 


0, 


j = 


1,.. 


,m, 


M r ^(h jy J(0)Ul(j)Ul(j^ 


-1)) 


y 


0, 


3 = 


1,.. 


,m, 


M r _^(^.y,/(0)U/(j)U/OH 


-1)) 


y 


0, 


3 = 


1,.. 


,m, 


MElli^j - !) 






0, 


3 = 


1,.. 


,m, 


MEjli - 1) 
Lj/Oy - U>ij) 






0, 


i = 


1,... 


,m, 






0, 


i,3 


= I,- 


. . , m 








1, 









with optimal value denoted inf Q r (and min Q r if the infimum is attained) . 



Then, with the notation 



Theorem 12. Let K C R n+m (compact) be the feasible domain of flOMRP A 
above: 

(a) inf Q r f px as i — > oo. 

(b) Let y r , be an optimal solution of the SDP relaxation (Q r ). If 

rankM r (y r ,/(0)) = rank M r ^ {y r , 1(0)) 
(8) rankM r (y r ,I(0)Ul(j)UJ(j + l)) = rank M,._, (y r , 1(0) U U I(j + 1)) j = I, m 

and i/rank(M r (y*,/(0) U U + 1)) n (J(j) U 7(j + 1)))) = 1 for all j then minQ r = p x . 

Moreover, let Aj := {(x* (j) , w* (j))} be the set of solutions obtained by the condition Q). Then, every 
(x*,w*) such that (x* , w*) i€ j^ = (x* (j) , w* (j)) for some Aj is an optimal solution of Problem MOMRF. 



Proof. The convergence of the semidefinite relaxation ( Q. r I was proved by Jibetean and De Klerk |T0j for 
a general rational function over a closed semialgebraic set. Here, we use that result applied to the rational 
function in Moreover, the index set of the indeterminates in the feasible set given by constraints 
(flJ-Q admits the decomposition I(k), k = 0. . . ,m that satisfies the running intersection property (see 
|17l (1.3)]) and therefore, the result follows by combining Theorem 3.2 in [17] and the results in [10]. □ 

The above theorem allows us to approximate and solve the original problem MOMRF up to any degree 
of accuracy by solving block diagonal (sparse) SDP programs which are convex programs for each fixed 
relaxation order r and that can be solved with available open source solvers as SeDuMi, SDPA, SDPT3 
[13], etc. 



4. Generalized Location Problems with rational objective functions 

This sections considers a wide family of continuous location problems that has attracted a lot of attention 
in the recent literature of location analysis but for which there are not common solution approaches. The 
challenge is to design a common resolution approach to solve them for different distances and in any finite 
dimension. 

We are given a set A = {ai, . . . ,a n } C K d endowed with an £ r -norm (here £ T stands for the norm 
i 

\\ x \\r = (Yli—i \ x i\ T ^j 7 f° r au x <= ^ rf ); an d a feasible domain K C R d , closed and semi-algebraic. The 

goal is to find a point x* € K C M. d minimizing some globalizing function of the distances to the set A. 
Here, we consider that the globalizing function is rather general and that it is given as a rational function. 
Some well-known examples are listed below (see e.g. [T], [3], 0) EH or [2S]) : 

n 

• f(u\, . . . , u n ) = \ui — Uj\, absolute deviation or envy problem. 

i<j 
n 

• f(ui, . . . , u n ) = V^(itj — u) 2 , variance problem. 

i=l 
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f(ui, . . . , u n ) = —L, obnoxious facility location. 



b 

• . . . , u„) = ^ — , Huff competitive location. 

J=1 1 + 

The main feature and what distinguishes location problem from other general purpose optimization prob- 
lems, is that the dependence of the decision variables is given throughout the norms to the demand points 
in A, i.e. \\x — aj|| r . In this section, we consider a generalized version of continuous single facility location 
problems with rational objective functions over closed semi-algebraic feasible sets. 

Let fj(u) := : E n h> 1, j = l,...,m be rational functions with qj(u) > for all j. Wc shall 

define the dependence of fj to the decision variable x E R d via u = (m, . . . ,u n ), whe re u, :R d ^ R, 
Ui{x) := \\x — di\\ Tl i — 1, . . . ,n. Therefore, the j-th component of the ordered median objective function 
of our problems reads as: 

fj(x) : R d i— > R 

x i v fj(x) :=/j-(||a;-oi|| T ,...,||a;-a n || T ). 

Consider the following problem: 

m 

(LOCOMRF) p x := min{^ Xj (x)f u) (x) : x e K }, 

j=i 

where: 

• K C R n satisfies Putinar's property, 

• t :— -, r, s G N, r > s and gcd(r, s) = 1. 

This problem does not reduce to the family MOMRF considered above since the dependence on the 
decision variable x is not given in the form of polynomials. Note that £ r -norms are not, in general, 
polynomials. 

To avoid this inconvenience, we introduce the following auxiliary problem. 



(9) Px = min ^2Xj(x)^2fi(u)wij 



= 1 i=l 
m 

s.t. ^ Wij = 1, for i = 1, . . . , m, 

3 = 1 



i=l 

m m 

X] w tjM u ) > W i3 + lM U )>3 = 1, • • • , "1, 

i=l i=l 

- w 4J = 0, for i, j = 1, . . . , to, 

(10) w|, > (x fe - a /d ) r , fe = 1, . . . , n, l = l,...,d, 

(11) v s kl > (a kl - xiY, k = 1, . . . , n, l = l,...,d, 

d 

u k = C^2 v ki) s , k = 1, . . . ,n, 
i=i 

m 

Wij £l, V i, j = 1, . . . ,to, 

Wfe; € R, fife € K, fc = 1, . . . , n, I = 1, . . . , d, 

x e K. 
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We note in passing that the above problem simplifies for those cases where r is even. In these cases, we 



can replace the two sets of constraints, namely (10) and (111 by the simplest constraint 

v ki = i x k - auf, Vfc, I. 

This reformulation reduces by (n x d) the number of constraints defining the feasible set. Moreover, these 
constraints do not induce semidefinite constraints in the moment approach but linear matrix inequalities 
which are easier to handle. Following the same scheme of the proof in Theorem [4] we get the following 
result, whose proof is left to the reader. 



Theorem 13. Let x be a feasible solution of ( LOCOMRF[ ) then there exists a solution (x, u, v, w) for 
such that their objective values are equal. Conversely, if (x,u,v,w) is a feasible solution for |$|) then there 
exists a solution (x) for (LOCOMRF) having the same objective value. In particular g\ =~Q\. Moreover, 
ifKc K d satisfies Putinar's property then K C R^" 1 + n ( d + 2 ) also satisfies Putinar's property. 

Now, we can prove a convergence result that allows us to solve, up to any degree of accuracy, the 
above class of problems. Let y = (y a ) be a real sequence indexed in the monomial basis (x@ u 1 v 5 of 
R[x,u,v,w] (with a = (£,7, S, Q G N d x N n x N nd x N™ 2 ). 

Let ho(x,u,v,w) := p\(x,u, v, w), and denote £j := |~(deg gj)/2~\ and Vj := [(degft 3 -)/2], where 



{gi, ■■■,££}, and {hi, 
K and K \ K in H. 

semidefinite programs 



(Qr) 



■ • ! "'3m+m 2 +2n((2+l) + l 

For 



} are, respectively, the polynomial constraints that define 



> r :— max{ max £ 



k=l, 



max i 

j=0,...,i+3m+m 2 +l 



introduce the hierarchy of 



mm 

y 
s.t. 



M r (y) 

M r _ 4fc (5fc,y) 



ho, 
to, 
to, 
= i, 



k = 1, 

i = i,. 



, 3m 



with optimal value denoted inf Q r (and minQ r if the iiifimuni is attained). 

Theorem 14. Let K C U d+m2+n( - d+2 '> (compact) be the feasible domain of Problem |s|). Let Q r be the 
semidefinite program |'Q r ). Then, with the notation above: 
(a) inf Q r f px as i — > oo. 



(b) Let y r be an optimal solution of the SDP relaxation Q r in |'Q r ). If 

rankM r (y r ) = rankM r _ ro (y r ) = t 

then minQ. r = P\ and one may extract t points {x*(k),U*(k),v*(k),w*(k)) k _ 1 C K, all global minimizers 
of the MOMRF problem. 



Proof. The convergence of the semidefinite relaxation Q r was proved by Jibetean and De Klerk [TU] for a 



general rational function over a closed semialgebraic set. Here, we apply this result applied to the rational 
function in ^ and therefore, the result follows. □ 

Here, we also observe that one can exploit the block diagonal structure of the problem since there is a 
sparsity pattern in the variables of formulation ([9|. The reader may note that the only monomials that 
appear in that formulation are of the form x a u? YYJLi v ii ^ or all * — L ■ ■ ■ j TO - Hence, a result similar to 



Theorem 14 also holds for the hierarchy (Q r ) of SDP applied to the location problem. Nevertheless, al- 



though we have used it in our computational test, we do not give specific details for the sake of presentation 
and because of the similarity with Theorem [14"} 

Example 15. We illustrate the above results with an instance of the well-known Weber problem with 
£ 3 -norm and for 20 random demand points in R 3 . Let A = {(0.0758,0.0540,0.5308), (0.7792,0.9340,0.1299), 
(0.5688, 0.4694, 0.0119), (0.3371, 0.1622, 0.7943), (0.3112, 0.5285, 0.1656), 

(0.6020, 0.2630, 0.6541), (0.6892, 0.7482, 0.4505), (0.0838, 0.2290, 0.9133), (0.1524, 0.8259, 0.5383), 
(0.9961, 0.0782, 0.4427), (0.1066, 0.9619, 0.0046), (0.7749, 0.8173, 0.8687), (0.0844, 0.3998, 0.2599), 
(0.8000, 0.4314, 0.9106), (0.1818, 0.2638, 0.1455), (0.1361, 0.8693, 0.5797), (0.5499, 0.1450, 0.8530), 
(0.5499, 0.1450, 0.8530), (0.4018, 0.0760, 0.2399), (0.1233, 0.1839, 0.2400)}. 
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Then, the problem consists of 



a£A 



S.t. 



x € 



The feasible region of the first SDP relaxation of this problem, which in this case is r = 2, contains 20 
moment matrices of size 36 x 36, 160 localizing matrices of size 8x8 and 36 equality constraints. The 
exact optimal solution is given by x = (0.426397,0.438730,0.455857) with optimal value f = 8.729976. 
We get with our approach, using SDPT3 13J, an optimal solution x* — (0.426397,0.438730,0.455857), 
for the first relaxation of the problem with optimal value f* = 8.729976. Thus, the relative error is 

e = i^yS = 2.199595 x 10~ 13 . 

For the same set of points, we consider a modification of the above problem by adding an extra nonconvex 
constraint: 



s.t. 



— 1x\ — 

x e 



x - a 3 



2x§ > 0, 



The exact optimal solution of this problem is x = (0.562304,0.266296,0.295262) with optimal value f = 
10.109333. The reader may note that the original solution x is not feasible for the new problem. Using 
our approach, again for the first relaxation order, we get x** — (0.562304, 0.266296, 0.295262) with optimal 



\r*-f\ = 
f 



5.801151 x 10" 



value /** = 10.109333. Hence, the relative error in this case is e 

We show in Figure [7] the feasible region of our problem as well as the demand points and the optimal 
solutions (the exact and the ones obtained with our relaxed formulations) of the problems. The demand 
points in A are represented by ' * ', the optimal solution, x* , of the SDP relaxation without the nonconvex 
constraint by ' m ' and the optimal solution, x** , of the SDP relaxation with the nonconvex constraint is 
depicted by ' • '. 




FIGURE 1 . Feasible region, demand points and optimal solutions of Example [l5j 
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In the following, we will apply this general methodology to get the reformulation of the most standard 
problems in Location Analysis (see Nickel and Puerto [35]) that will be later the basis of our computa- 
tional experiments: minisum (Weber) and minimax (center), fc-centrum, (fci, fc2)-trimmed mean and range 
problems. 



4.1. Weber or median problem. In the standard version of the Weber problem, we are given a set of 
demand points {ai, . . . , a n } in M. d and a set of non-negative weights wi,...,w„ and one looks for a point 
x* minimizing the weighted Euclidean distance from the demand point. In other words, the problem is: 



min } LUi 

x£R d ' 
i—1 



This problem has been largely studied in the literature of Location Analysis and perhaps its most well- 
known algorithm is the so called Weiszfeld algorithm (see [38] ) . This problem is a convex one and Weiszfcld 
algorithm is a gradient type iterative algorithmic scheme for which several convergence results are known. 



Here, we observe that this problem corresponds to a very particular choice of the elements in ( LOCOMRF I : 
A = (1, . . . , 1), fi(u) — u>iU and r = 2, s = 1. Furthermore, the general formulation (LOCOMRF I simpli- 
fies since there is no actual sorting. Therefore, we can avoid many of our instrumental variables, namely, 
the problem can be cast into the form: 



n 

mm LU.iZ, 
1=1 

d 

S 

3=1 



i=l 



(WP) £>| + z? <M,i = l,...,7i, 

Zi > 0, i = 1, ... ,n, 
x e R d . 



4.2. The minimax or center problem. The minimax location problem looks for the location of a server 
x G M. d that minimizes the maximum weighted distance to a given set of demands points {ai, . . . , a„} in 
M. d . Formally, the problem can be stated as: 



min max u)i\\x — a^\%, 

kGR c! i=l,...,n 



for some weights w\, . . . , u n > 0. 

Once more, this problem has been extensively analyzed in the literature of Location Analysis and 
the most well-known algorithms to solve it are those by Elzinga-Hearn (only valid in M 2 with Euclidean 
distance) and Dyer [51 [7] and Megiddo [23] which are polynomial in fixed dimension. Again, we observe that 
this problem corresponds to a very particular choice of the elements in (LOCOMRF): A = (1,0, ... , 0), 
fi(u) = u)iU and r = 2, s = 1. In this case, the general formulation (LOCOMRF ) simplifies and therefore, 
we can avoid many of our instrumental variables, namely, the problem can be formulated as: 
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min t 

d 

s.t. zf = ^2(xj - a^) 2 , i = l,...,n, 
3=1 

LOiZi <t, i = l,...,n, 

d 

(CP) ^x] + z 2 +t 2 < M, i = 1,.. 

3=1 

*, Zi,>0, i=l,...,n, 
xeK. 

4.3. The k-centrum problem. The fc-centrum location problem consists of finding the point x* that 
minimizes the sum of the k largest distances with respect to a given set of demands points {a%, . . . , a n } in 
R d . Formally, the problem can be stated as: 

min max da)(x), 

where d^(x) — \\x — dcr(t)|| for a permutation a such that d a (i){x) > . . . > d a ^(x). This problem has been 
considered in several papers and textbooks (see [25], [I]). Currently, there exist few approaches to solve it 
in the plane (i.e. d = 2) and with the Euclidean norm that do not extend further to higher dimension nor 
other norms (see [HGIJIM]). The objective function of this problem is described by a vector of A-parameters 



A = (1, . . . , 1, 0, . . . , 0), fi(u) = u, r = 2, s = 1. Using the result in the reformulation (kC I the problem 
can be restated as: 



g := min kt + r,; 

i=i 
d 

s - t - z i = - a ij) 2 , * = i, ... ,n, 

3=1 

(kCP) t + rj > i = 1, . . . , n, 

^;c 2 + z 2 +r 2 < M, i = l,...,n, 

3=1 

i,r,,Zi > 0, i = 1, ... ,ri, 
jc G R d . 

4.4. The (fci, ^J-trimmed-mean problem. The (hi, fc2)-trimmed-mean location problem looks for a 
point x* that minimizes the sum of the central distances, once we have excluded the fc 2 closest and the k\ 
furthest. Formally, the problem is: 

min > df,-i(x), 

i=k 1 + l 

where di^{x) — \\x — H2 f° r a permutation cr such that d a ^(x) > ... > d a / n \(x). This problem has 
been considered in several papers and textbooks (see [5S], [1]). Currently, there exists two approaches 
to solve it in the plane (i.e. d = 2) and with the Euclidean norm that do not extend further to higher 
dimension nor other norms (see 016!). The objective function of this problem, in terms of the elements in 



(LOCOMRF), is described by a vector of A-parameters A = (0, . . . , 0, 1, . . . , 1, 0, . . . , 0), fi(u) = u, r = 2, 



s = 1. Here, we could apply the general formulation derived from (LOCOMRF). Nevertheless, that 



approach needs many decision variables which affects the sizes of the problems to be handled. Rather than 
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the general formulation, we present here an alternative problem, based on (kTr), which takes advantage of 



the particular structure of this problem and reduces the number of variables needed for its representation. 
We consider the problem: 



(TMP) 



min (n - k 2 )t + ^ r t - ^ u i z l 

i=l i=l 

d 

s.t. zf = ^(xj - fl tJ ) 2 , i = l,...,'i 



3=1 



^ Ui = ki , 

i=i 

m(ui - 1) = 0, i = 1, . . . ,n, 

t + r,; > Zi . 1=1 



! " ' 1 '"J 



^ + z 2 + t 2 + w 2 + r 2 < M, i = 1 

3 = 1 



) * ' ' 5 5 



a; G 



4.5. The range problem. The last problem that we address in our computational experiments is the 
range location problem. This problem consists of minimizing the difference (range) between the maximum 
and minimum distances with respect to a given set of demands points {eti, . . . ,a n } in K d (see El I2S"] ). 
Formally, the problem can be stated as: 



mm max \\x — aj|| 2 — mm \\x — dj||2 

a;gB d i=l,...,n t=l,...,n 



This problem corresponds to the following choice of the elements in (LOCOMRF I: A = (1, 0, . . . , 0, —1), f%(u) 
u and r = 2, s = 1. A simplified reformulation of the problem reduces to: 



min z — t 

d 

s - 1 - z l = ^2( x j - a y) 2 ) « = 1, - ■ • ,n, 
J=l 

(RP) i < Zi < z, i = 1, . . . , n, 

d 

x 2 + z 2 + t 2 + z 2 < M, i = l,...,n, 

3 = 1 

t,z,Zi>0, i = l,...,n, 
x G M d . 



5. Computational Experiments 

A series of computational experiments have been performed in order to evaluate the behavior of the 
proposed methodology. Programs have been coded in MATLAB R2010b and executed in a PC with an 
Intel Core i7 processor at 2x 2.93 GHz and 8 GB of RAM. The semidefinite programs have been solved by 
calling SDPT3 4.0Q3]. 

We run the algorithm for several well-known continuous location problems: Weber problem, center 
problem, k-center problem, trimmed-mean problem and range problem. For each of them, we obtain the 
CPU times for computing solutions as well as the gap with respect to the optimal solution obtained with 
the battery of functions in optimset of MATLAB or the implementation by [6]. 
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With regard to computing the accuracy of an obtained solution, we use the following measure for the 
error (see |37|): 

| the optimal value of the SDP — fopt| 
( } £obj = max{l,fopt} ' 

where fopt is the optimal objective value for the problem obtained with the functions in optimset or the 
implementation by [5] • 

We have organized our computational experiments in five different problems types that coincide with 
those described previously in sections 4.1|4.5 Our test problems are generated to be comparable with 



previous results of some algorithms in the plane but, in addition, we also consider problems in R . Thus, 
we report on randomly generated points on the unit square and in the unit cube. Depending on the 
problem, we have been able to solve different problem sizes. In all problems, we could solve instances with 
at least 500 points for planar and 3-dimensional problems and with an average accuracy higher than 10 -5 . 
(We remark that for instance we could solve instances of more than 1000 points for Weber and center 
problems with high precisions.) 

Our goal is to present the results organized per problem type, framework space (M 2 or ]R 3 ) and relaxation 
order. We report for Weber problem on the first two relaxations to show that raising relaxation order one 
gains some extra precision (as expected) at the price of higher CPU times. In spite of that, the considered 
problems seems to be very well-approximated even with the first relaxation (as shown by our results). For 
this reason, we only report results for relaxation order r — 2 for the remaining problem types, namely 
center, fc-centrum, range and trim-mean. 

The results in our tables, for each size and problem type, are the average of ten runs. In all cases our 
tables are organized in the same way. Rows give the results for the different number of demand points 
considered in the problems. Column n stands for the number of points considered in the problem, CPU time 



is the average running time needed to solve each of the instances, e a bj gives the error measure (see 12 1. 
The final block of 3 columns informs on the sizes of the SDP problems to be solved: #Cols, #Rows and 
°/ NonZero represent, respectively, the number of columns, rows and the percentage of nonzero entries of 
the constraint matrices of the problems to be considered. 

We tested problems with up to 500 demands points (except for Weber problem where we considered 
1000 demands points) randomly generated in the unit square and the unit cube. We move n between 10 and 
500 (or 1000 for Weber problem) and ten instances were generated for each value of n. The first relaxation 
of the problems was solved in all cases. For the k-centrum problem type we considered three different k 
values to test the difficulty of problems with respect to that parameter, k — [O.ln], [~0.5n~|, |~0.9n] (tables 
|4]and[5|. 

Tables [T][7] show the averages CPU times and gaps obtained. Table [T] summarizes the results of the 
Weber problems. We remark that problems with up to 1000 demand points on the plane are solved with 
the first relaxation in few seconds and with accuracy higher than 10~ 4 . Raising the relaxation order, we 
improve accuracy till 10 -6 at the cost of multiplying CPU time by a factor of 8. Table [2] refers to Weber 
problem in the 3d space. Results are similar although precision is higher when considering the second 
relaxation order. Table [3] reports the results for the center problem on the plane and the 3d-space. CPU 
times are slightly larger than for the Weber problem but accuracy are also better specially for sizes up to 
100 demand points. Tables [4] and [5] are devoted to show the behavior of our approach for three different 
k values of the fc-centrum problem (Table [4] in K 2 and Table [5] in K 3 ). We observe that for small values 
of k, i.e. k = [O.ln] or [~0.5n] the fc-centrum is slightly harder than for values closer to n. The remaining 
factors behave similarly to those in Weber or center problems. Table [6] reports the results for the range 
problem. The behavior of these problems is similar to that of the fccentrum problems both in CPU time 
and accuracy. Finally, Table [7] summarizes the results for the trimmed- mean problems. These are the 
harder problems among the five considered problem types. We are able to solve similar sizes with similar 
accuracies using the first order relaxation. However, CPU times are significantly higher than for the other 
problem types. These results show that this methodology can be efficiently applied to solve medium to 
large sized location problems. 

^From our tables we conclude that Weber problem is the simplest one whereas the trimmed-mean 
problem is the hardest one, as expected. We remark that CPU times increase linearly with the number 
of points in all problem types. A linear regression between these times and the number of points gives 
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a regression coefficient i?-squared (coefficient of determination of the regression) greater than 0.98 for all 
the problems. Therefore, this shows a linear dependence, up to the tested sizes, between problem sizes 
and CPU times for solving the corresponding relaxations. Observe that the sizes of the matrices in the 
SDP relaxations increase exponentially with the number of points. Nevertheless, the percentage of nonzero 
elements in the constraint matrices decreases very slowly (hyperbolically) when increasing the size (number 
of points) of the problems. 

6. Conclusions 

We develop a unified tool for minimizing weighted ordered averaging of rational functions. This approach 
goes beyond a trivial adaptation of the general theory of moments-sos since ordered weighted averages of 
rational functions are not, in general, neither rational functions nor the supremum of rational functions so 
that current results cannot directly be applied to handle these problems. As an important application we 
cast a general class of continuous location problems within the minimization of OWA rational functions. 
We report computational results that show the powerfulness of this methodology to solve medium to large 
continuous location problems. 

This new approach solves a broad class of convex and non convex continuous location problems that, up 
to date, were only partially solved in the specialized literature. We have tested this methodology with some 
medium to large size standard ordered median location problems in the plane and in the 3-dimcnsional 
space. Our goal was not to compete with previous algorithms since most of them are either problem 
specific or only applicable for planar problems. However, in all cases we obtained reasonable CPU times 
and high accuracy results even with first relaxation order. Our good results heavily rely on the fact that 
we have detected sparsity patterns in these problems reducing considerably the sizes of the SDP object to 
be considered. 

The two main lines for further research in this area would be to increase both the sizes and the classes 
of problems efficiently solved. These goals may be achieved by improving the efficiency of available SDP 
solvers and/or by finding alternative formulations that take advantage of new sparsity and symmetry 
patterns. 
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Table 1. Computational results for planar Weber problem and first and second relaxation. 
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Table 2. Computational results for Weber problem in E. 3 and first and second relaxation. 
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Table 3. Computational results for center problem in M 2 and M 3 and first relaxations. 
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0.036% 


60 


25.77 


0.00000077 


25.41 


0.00000283 


20.21 


0.00000806 


39420 


5494 


0.030% 


70 


28.01 


0.00000204 


31.02 


0.00000234 


25.07 


0.00000192 


45990 


6404 


0.026% 


80 


37.25 


0.00000085 


31.48 


0.00000044 


30.66 


0.00000220 


52560 


7314 


0.023% 


90 


47.16 


0.00000062 


41.07 


0.00000765 


33.92 


0.00000086 


59130 


8224 


0.020% 


100 


53.68 


0.00000084 


41.42 


0.00000065 


39.49 


0.00000188 


65700 


9134 


0.018% 


150 


86.48 


0.00000089 


68.48 


0.00000056 


65.95 


0.00000059 


98550 


13684 


0.012% 


200 


123.02 


0.00000056 


96.40 


0.00000075 


88.10 


0.00000275 


131400 


18234 


0.009% 


250 


149.26 


0.00003681 


135.67 


0.00000071 


113.68 


0.00000161 


164250 


22784 


0.007% 


300 


180.38 


0.00000408 


161.84 


0.00000081 


146.22 


0.00000349 


197100 


27334 


0.006% 


350 


223.27 


0.00003013 


193.31 


0.00003623 


176.46 


0.00000151 


229950 


31884 


0.005% 


400 


260.27 


0.00000079 


225.07 


0.00003689 


201.01 


0.00000376 


262800 


36434 


0.005% 


450 


290.23 


0.00004512 


272.55 


0.00000097 


237.23 


0.00000168 


295650 


40984 


0.004% 


500 


345.93 


0.00000224 


310.19 


0.00000119 


269.99 


0.00000200 


328500 


45534 


0.004% 



Table 4. Computational results for planar fc-ccntrum problems and first relaxation (r = 2). 
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k = 


[0.1 n] 


k = 


[0.5 n] 


k = 


[0.9 n] 


Sizes 


n 


CPU time 


£obj 


CPU time e Q bj 


CPU time e ob j 


#Cols 


#Rows 


7oNonZero 


10 


7.06 


0.00041340 


5.85 


0.00000039 


6.05 


0.00000168 


10780 


1469 


0.114% 


20 


16.40 


0.00000950 


15.42 


0.00000095 


16.30 


0.00000019 


21560 


2869 


0.059% 


30 


27.63 


0.00001682 


23.72 


0.00000028 


27.12 


0.00000132 


32340 


4269 


0.039% 


40 


45.25 


0.00000075 


42.31 


0.00000086 


37.38 


0.00000077 


43120 


5669 


0.030% 


50 


54.39 


0.00000282 


53.66 


0.00000026 


51.94 


0.00000087 


53900 


7069 


0.024% 


60 


63.16 


0.00000259 


59.34 


0.00000091 


63.91 


0.00000065 


64680 


8469 


0.020% 


70 


85.17 


0.00000144 


81.32 


0.00000258 


74.24 


0.00000079 


75460 


9869 


0.017% 


80 


106.65 


0.00000326 


83.96 


0.00000044 


88.76 


0.00000158 


86240 


11269 


0.015% 


90 


114.38 


0.00000209 


93.85 


0.00000100 


103.56 


0.00000092 


97020 


12669 


0.013% 


100 


122.01 


0.00000088 


109.17 


0.00000224 


118.03 


0.00000067 


107800 


14069 


0.012% 


150 


235.10 


0.00000073 


211.54 


0.00000890 


187.51 


0.00000135 


161700 


21069 


0.008% 


200 


305.51 


0.00002407 


255.54 


0.00007106 


284.80 


0.00000157 


215600 


28069 


0.006% 


250 


403.89 


0.00000519 


348.32 


0.00004300 


357.79 


0.00000143 


269500 


35069 


0.005% 


300 


492.04 


0.00046130 


433.69 


0.00007630 


471.78 


0.00000174 


323400 


42069 


0.004% 


350 


529.61 


0.00041229 


484.87 


0.00000058 


448.60 


0.00001791 


377300 


49069 


0.003% 


400 


619.97 


0.00000091 


585.93 


0.00000055 


523.81 


0.00000829 


431200 


56069 


0.003% 


450 


705.99 


0.00048727 


693.77 


0.00000037 


580.06 


0.00004327 


485100 


63069 


0.003% 


500 


817.75 


0.00012138 


789.77 


0.00000087 


664.94 


0.00000318 


539000 


70069 


0.002% 



Table 5. Computational results for fc-centrum problems in R 3 and first relaxation (r = 2). 







M 3 


n 


CPU time e ob j 


#Cols 


#Rows 


7oNonZero 


CPU time 


e Q bj 


#Cols 


#Rows 


°/ NonZero 


10 


2.96 


0.00007519 


6060 


629 


0.252% 


5.68 


0.00001997 


10080 


965 


0.164% 


20 


7.04 


0.00001750 


12120 


1189 


0.133% 


18.45 


0.00015758 


20160 


1805 


0.088% 


30 


13.94 


0.00098322 


18180 


1749 


0.091% 


35.37 


0.00028187 


30240 


2645 


0.060% 


40 


14.53 


0.00002124 


24240 


2309 


0.069% 


35.77 


0.00032049 


40320 


3485 


0.045% 


50 


24.49 


0.00004314 


30300 


2869 


0.055% 


65.80 


0.00051293 


50400 


4325 


0.037% 


60 


23.49 


0.00047832 


36360 


3429 


0.046% 


59.19 


0.00005082 


60480 


5165 


0.031% 


70 


34.87 


0.00003903 


42420 


3989 


0.040% 


68.46 


0.00006841 


70560 


6005 


0.026% 


80 


38.69 


0.00026693 


48480 


4549 


0.035% 


79.54 


0.00003016 


80640 


6845 


0.023% 


90 


42.34 


0.00042121 


54540 


5109 


0.031% 


90.76 


0.00017468 


90720 


7685 


0.021% 


100 


58.36 


0.00052427 


60600 


5669 


0.028% 


97.26 


0.00015535 


100800 


8525 


0.019% 


150 


65.04 


0.00021457 


90900 


8469 


0.019% 


159.41 


0.00094711 


151200 


12725 


0.012% 


200 


98.23 


0.00041499 


121200 


11269 


0.014% 


197.66 


0.00040517 


201600 


16925 


0.009% 


250 


131.42 


0.00033959 


151500 


14069 


0.011% 


274.14 


0.00057559 


252000 


21125 


0.007% 


300 


159.87 


0.00014556 


181800 


16869 


0.009% 


322.21 


0.00036845 


302400 


25325 


0.006% 


350 


169.29 


0.00003661 


212100 


19669 


0.008% 


393.80 


0.00096204 


352800 


29525 


0.005% 


400 


167.74 


0.00123896 


242400 


22469 


0.007% 


361.12 


0.00022448 


403200 


33725 


0.005% 


450 


218.70 


0.00207328 


272700 


25269 


0.006% 


513.55 


0.00044016 


453600 


37925 


0.004% 


500 


228.68 


0.00438388 


303000 


28069 


0.006% 


554.94 


0.00028013 


504000 


42125 


0.004% 



Table 6. Computational results for range problem in M 2 and K 3 and first relaxation. 
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n 


CPU time 


e bj 


#Cols 


#Rows 


%NonZero 


CPU time 


£ obj 


#Cols 


#Rows 


%NonZero 


10 


5.31 


0.00017041 


11760 


1784 


0.087% 


14.09 


0.00000197 


18080 


2669 


0.059% 


20 


12.39 


0.00000619 


23520 


3534 


0.044% 


33.85 


0.00047792 


36160 


5269 


0.030% 


30 


18.11 


0.00020027 


35280 


5284 


0.029% 


49.16 


0.00000670 


54240 


7869 


0.020% 


40 


30.39 


0.00035248 


47040 


7034 


0.022% 


73.13 


0.00001450 


72320 


10469 


0.015% 


50 


36.04 


0.00181487 


58800 


8784 


0.018% 


98.17 


0.00001624 


90400 


13069 


0.012% 


60 


49.16 


0.00085810 


70560 


10534 


0.015% 


131.38 


0.00003143 


108480 


15669 


0.010% 


70 


60.57 


0.00012995 


82320 


12284 


0.013% 


161.25 


0.00004420 


126560 


18269 


0.009% 


80 


73.54 


0.00092073 


94080 


14034 


0.011% 


188.51 


0.00012265 


144640 


20869 


0.008% 


90 


76.12 


0.00040564 


105840 


15784 


0.010% 


203.06 


0.00011847 


162720 


23469 


0.007% 


100 


91.26 


0.00218668 


117600 


17534 


0.009% 


220.68 


0.00011032 


180800 


26069 


0.006% 


150 


153.31 


0.00814047 


176400 


26284 


0.006% 


400.37 


0.00026203 


271200 


39069 


0.004% 


200 


257.23 


0.00032380 


235200 


35034 


0.004% 


552.19 


0.00056138 


361600 


52069 


0.003% 


250 


339.72 


0.00051519 


294000 


43784 


0.004% 


659.01 


0.00046219 


452000 


65069 


0.002% 


300 


326.52 


0.00225994 


352800 


52534 


0.003% 


884.40 


0.00038481 


542400 


78069 


0.002% 


350 


410.32 


0.00047898 


411600 


61284 


0.003% 


955.53 


0.00061467 


632800 


91069 


0.002% 


400 


582.36 


0.00047130 


470400 


70034 


0.002% 


1165.79 


0.00058261 


723200 


104069 


0.002% 


450 


631.58 


0.00060180 


529200 


78784 


0.002% 


1931.76 


0.00081711 


813600 


117069 


0.001% 


500 


685.79 


0.00079679 


588000 


87534 


0.002% 


9151.90 


0.00063861 


904000 


130069 


0.001% 



Table 7. Computational results for trimmed mean problem with fcl = k2 = [0.20 n] in W 2 and M 3 and first relaxation. 
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